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Abstract. This article is devoted to incompressible Euler equations (or to Navier- 
Stokes equations in the vanishing viscosity limit). It describes the propagation of 
quasi- singularities. The underlying phenomena are consistent with the notion of a 
cascade of energy. 

Resume. Get article etudie les equations d'Euler incompressible (ou de Navier- 
Stokes en presence de viscosite evanescente) . On y decrit la propagation de quasi- 
singularites. Les phcnomenes sous-jacents confirmcnt I'idee selon laquelle il se pro- 
duit une cascade d'energie. 

1 Introduction. 

Consider incompressible fluid equations 

{£) Stu + (u • V) u + Vp = , div u = , {t, x) £ [0, T] x M*^ , 

where u = *(u^,--- , u*^) is the fluid velocity and p is the pressure. The 
structure of weak solutions of {£) in d— space dimensions with d > 2 is a 
problem of wide current interest [4]- [23]. The questions are how to describe 
the phenomena with adequate models and how to visualize the results in 
spite of their complexity. We will achieve a small step in these two directions. 

According to the physical intuition, the appearance of singularities is 
linked with the increase of the vorticity. Along this line, we have to mark 
the contributions [2] and [9]. Interesting objects are solutions which do not 
blow up in finite time but whose associated vorticities increase arbitrarily 
fast. These are quasi- singularities. Their study is of practical importance. 

Typical examples of quasi-singularities are oscillations. This is a well- 
known fact going back to [3]- [24]. The works [3] and [24] rely on phenomeno- 
logical considerations and engineering experiments. Further developments 
are related to homogenization [12]-[13], compensated compactness [11]-[16] 
and non linear geometric optics [6]- [7]- [8]. 
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In [11], DiPerna and Majda show the persistence of oscillations in three 
dimensional Euler equations. To this end, they select a parameter s G ]0, 1] 
and look at 

(1.1) ul{t, x) := *(g(a:2, e"^ X2), 0, h(xi - g(a:2, e"^ X2) t, X2, ^2)) 

where g{x2,6) and h(xi,X2, 9) are smooth bounded functions with period 1 
in 9. They remark that the functions are exact smooth solutions of {£) 
and they let e goes to zero. Yet, this construction is of a very special form. 
First, it comes from shear layers (steady 2-D solutions) as 

u^(t,x) = u|(0,x) = *(g(x2,e-ix2),0) G . 

Secondly, it involves a phase (pQ{t,x) = X2 which does not depend on e. 
Of course, this is a common fact [10]- [18]- [19]- [26] when dealing with large 
amplitude high frequency waves. Nevertheless, this is far from giving a 
complete idea of what can happen. 

Our aim in this paper is to develop a theory which allows to remove the 

two restrictions mentioned above. Fix b = {l,N) G where the integers I 
and N are such that < I < N . Introduce the geometrical phase 

/I ^ 

iflit.x) := ipoit,x) + J2k=i "Pkit.x). 

In the section 2, we state the Theorem 2.1 which provides with approximate 
solutions defined on the interval [0, T] with T > and having the form 

(12) "'^^^'''^ = ,0(*,^) 

= uo(t,x) + Ef=i £^ Uk{t,x,e-^^l{t,x)) 

where the smooth profiles 

Uk{t, X, 9) = \Ul Ui){t, X, 0) G 1 < A; < AT , 

are periodic functions of G M/Z. We assume that 

3{t,x,e) G [0,r] X M'' X T; deUi{t,x,0) 7^0. 

We say that the family {u^je is a weak, a strong or a turbulent oscillation 
according as we have respectively Z = 1, Z = 2 or Z > 3. 

The order of magnitude of the energy of the oscillations is e~. Com- 
pute the vorticities associated with the functions u^. These are the skew- 
symmetric matrices $7^ = {^lj)i<i,j<d where 

J7g(t,x) := {djul'-d,ul^){t,x) 

= Ek=i e^"' (^M deUl-d^^l deUl){t,x,e-^cplit,x)) 

+ {djul,-diui){t,x) + ST {djUl-diUi){t,x,e-^^l{t,x)). 
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The principal term in ft^ is of size s'i~^. When Z > 2, no uniform estimates 
are available on the family {^^^jegjo,!]- In particular, if d = 3, there is no 
uniform control on the enstrophy 

/o^/iR3 K(i,a;)|2 dtdx, io'^{t,x) := (VAup(t,x) = n'^{t,x). 

We see here that strong and turbulent oscillations are examples of quasi- 
singularities. Observe that the expansion (1.2) involves a more complicated 
structure than in (1.1) though the corresponding regime is less singular. 

The BKW analysis reveals that the phase shift ipi and the terms ipk 
with 2 < A; < Z — 1 play different parts. The role of ipi is partly examined 
in the articles [6] and [7] which deal with the case / = 2. When / > 3, the 
phenomenon to emphasize is the creation of the (pk with 2 < k < I — 1. 
Indeed, suppose that 

^2(0,0= ••• =^i-i(0, = 0, l>3. 

Then, generically, we find 

3te]0,T]; ^2{t,-)^0, <^,_i(t,-)#0. 

Now starting with large amplitude waves (this corresponds to the limit case 

I = +oo) that is 

the description of u^(t, •) on the interval [0, T] with T > needs the intro- 
duction of an infinite cascade of phases ipk- The scenario is the following. 
Oscillations of the velocity develop spontaneously in all the intermediate fre- 

k 1 

quencies e'i~ and in all the directions V(pk{t,x). This expresses turbulent 
features in the flow. 

The family {u^Isgjo,!] is e— stratifled [19] with respect to the phase 
with in general ^ ipQ. The presence in cpg of the non trivial functions 
ifk is necessary and sufficient to encompass the geometrical features of the 
propagation. It has various consequences which are detailed in the section 
3. It brings informations about microstructures, compensated compactness 
and non linear geometric optics. It also confirms observations made in the 
statistical approach of turbulences [14]- [22]. 

The chapter 4 is devoted to the demonstration of Theorem 2.1. Because 
of closure problems, the use of the geometrical phase Pg does not suffice 
to perform the BKW analysis. Among other things, adjusting phases ipk 
with I < k < N must be incorporated in order to put the system of formal 
equations in a triangular form. 
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The expressions are not exact solutions of Euler equations, yielding 
small error terms as source terms. The matter is to know if there exists 
exact solutions which coincide with u^(0, •) at time t = 0, which are defined 
on [0, T] with T > 0, and which are close to the approximate divergence free 
solutions u^. This is a problem of stability. 

The construction of exact solutions requires a good understanding of the 
different mechanisms of amplifications which occur. In the subsection 5.1, 
we make a distinction between obvious and hidden instabilities. 

The obvious instabilities can be detected by looking at the BKW analysis 
presented before. They imply the non linear instability of Euler equations 
(Proposition 5.1). They need to be absorbed a dependent change of variables 
which induces a defect of hyperbolicity. The hidden instabilities can be 
revealed by soliciting this lack of hyperbolicity. They require to be controled 
the addition of dissipation terms. 

In the subsection 5.2, we look at incompressible fiuids with anisotropic vis- 
cosity. This is the framework of [5] though we adopt a different point of 
view. We consider strong oscillations. We show (Theorem 5.1) that exact 
solutions corresponding to 11^2 n) ^^i^^ on some interval [0, T] with T > 
independent on e G ]0, 1]. 
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2 Euler equations in the variables {t,x). 



The description of incompressible flows in turbulent regime is a delicate 
question. No systematic analysis is yet available. However, special appro- 
ximate solutions with rapidly varying structure in space and time can be 
exhibited. Their construction is summarized in this chapter 2. 



2.1 Notations. 

• Variables. Let T G M+. The time variable is t e [0,T]. Let d G N \ {0, 1}. 
The space variables are (x, (9) G M'' X T where T := R/Z. Mark the bah 

B{0, R] := { X G M"' ; jxp := J^Li xf < R] , G M+ . 

The state variables are the velocity field u = ,u'^) G M"^ and the 

pressure p G M. Given {u,u) G (M*^)^, define 

u ■ u := "U* -u* , := u ■ u , u := {u^ u^)i<i,j<d ■ 

The symbol S'^ is for the set of positive definite quadratic form on W^. An 
element q e Sf. can be represented by some dx d matrix {c{ij)i<i,j<d- 

• Functional spaces. Distinguish the expressions u(t, x) which do not depend 
on the variable 6 from the expressions u{t,x,6) which depend on 6. The 
boldfaced type u is used in the first case whereas the letter u is employed 
in the second situation. 



Note C^°°([0,r] X W) the space of functions in [0,r] x R'^ with bounded 
continuous derivatives of any order. Let m G N. The Sobolev space H"^ is 
the set of functions 

such that 

II u Wjjm := E.ez (1 + \k\'r /m. (1 + l^T < ^ 

where 

J^(u)(0 = n{0 ■■= (27r)-t J^a e-'-< u(x) dx , ^ G M'^ . 
With these conventions, the condition u G H'^ means simply that 

\\^\?Hr.:= !^a{l + \er |u(6Ne<oo. 
Define 

HJp := {u; din G L2([0, T]; /f'"-^) , V j G {0, • • • , m} } , 
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WJp := [u-ue CJ([0, T]-H^-i) , Vj e {0, ■ ■ ■ , m} } , 
with the corresponding norms 

II W'ujp ■~ ^j=0 fo II ll/f™ '^^ ^ 

II u llwj? := suptg[o,T] EjLo II ^"(^' •) Hi?"" ■ 
Consider also 

-"oo •— 1 lTeR+ -"T 1 -"T 1 ImeN -"T ' -"oo •— I lTeK+ -"T ' 

:= nTeM+ >V?^ ' ■■= HmeN > := nTeM+ '^T • 

When m = 0, replace with L^. Any function u e can be decomposed 
according to 

u{t,x,6) = {u){t,x) + u*{t,x,9) = u{t,x) + u*{t,x,6) 
where 

{u){t,x) = u{t,x) := Jrj, u{t,x,9) d9 . 

Let r be the symbol of any of the spaces H"^, Hip, W™, • • • defined before. 
In order to specify the functions with mean value zero, introduce 

r* := {ueV; u = 0}. 

Mark also 

suppj. u* := closure of { x G M'' ; || u*{x, ■) ||l2(t)7^ } . 

• Differential operators. Note 

dt = do := d/d t, dg = da+i := d/d 6 , 

dj --d/dxj, j £{!,■■■ ,d}, 

V:={di,--- ,dd), A:=A, + dl = d^ + --- + d^ + d',. 

Let u e W^. Define 

u-V := di-\ \-u'^ dd, 

div u := diu^ + • • • + d^u'^ , 

div (n (g) u) := Yl'i=i ^{djiu^ u^),-- - , dj{u^ u'^)) G R"^ . 

Employ the bracket < •, • >h for the scalar product in the Hilbert space 
H. Note C{E;F) the space of linear continuous applications T : E — > F 
where E and F are Banach spaces. The symbol C(E) is simply for C{E; E). 
Introduce the commutator 

[A-B] := AoB-BoA, {A,B) e C{Ef . 
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Let r G Z. The operator T is in if 

II \\c{Hjp+^-H;p) ^ °° ' ^ ^ • 

Let £o > 0. The family of operators {T^}^ G is in if 

SUPee]0,£o] II ll£(i/-+'-;Jf-) < °° ' Vm € N . 
Consider a family {f'}e € (>V°°)]0'^<'1. We say that = 0{e'') if 

suPee]o,£o] II llw?? < oo , Vm G N. 

Given a family {f^}^ G we say that {Fj^ = O(e0 if 

sup,g]o,eo] ||fMlw-<oo, VmeN. 

Observe that the two preceding definitions have very different significations 
according as we use the letter / or the boldfaced type f . In particular, the 
second inequalities correspond to e— stratified estimates. The families {f^}e 
or {f^je are C){e°°) if they are C){s^) for all r G M. 

2.2 Divergence free approximate solutions in {t,x). 
• A first result. Select smooth functions 

uoo €H^, ^00 e (M<^) , V^oo e (K'^) . 
Suppose that 

3c>0; \V(poo{x)\ > 2 c, VxGM'^'. 
For T > small enough, the equation {£) associated with 

uo(0,x) = uoo(a;), VxgM'^ 
has a smooth solution uo(i,x) G W|?. Solve the eiconal equation 
(ei) dt^o + (uo • V) (^0 = , {t, x) e [0, T] x M'^ 
with the initial data 

If necessary, restrict the time T in order to have 

(2.1) \Vipo{t, x)\ > c, V (t, x) G [0, T] X M"^ . 

Call Ilo{t,x) the orthogonal projector from onto the hyperplane 
V(po(t,x)-^ := (ugM''; u ■ Vipo{t,x) = O} . 
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Theorem 2.1. Select any b = {l,N) G such that < Z(3 + ^) < iV. 

Consider the following initial data 

C/,*o(x, 9) = Uo{0, x) U*^^{x, e)eH^, l<k<N, 
Uko{x)eH°°, l<k<N, 
(Pkoix) eH°° , l<k<l-l. 

First, there are finite sequences {C^fe}i<fe<Ar o,nd {-Pjt}i<fe<Ar with 
Uk{t,x,9) eW^ , Pk{t,x,e) eW^ , l<k<N, 

and a finite sequence {'Pk}i<k<l-i with 
iPk{t,x)eW^, l<k<l-l, 

which are such that 

no(o,x)c/*(o,x,0) = no(o,x)c/*o(x,0), i< A;< at, 

Uk{0,x) = Uko{x), l<k<N, 

cpk{0,x) = ^pko{x) , l<k<l-l. 
Secondly, there is £o G ]0, 1] and correctors 

cu^{t, x)eW^, cpf (t, x)eW^, ee ]0, eo] , 
which give rise to families satisfying 

Then, all these expressions are adjusted so that the functions and 
defined according to 

^2 2) := uo{t,x) + Ylk=i Uk{t,x,e-^ ^pl{t,x)) +cul{t,x) 

pl{t,x) ■=Po{t,x) + J2k=i Pk{t,x,e^^ ipl{t,x)) +cpl{t,x) 

where (pg{t,x) is the geometrical phase 

(2.3) ipl{t,x) := ipQ{t,x) + Ei=\ Vk{t,x) 

are approximate solutions of {£) on the interval [0, T] . More precisely 

dtMl + (4 • V)u^ + Vpl = , div = , = 0{eT-^-i) . 
• Some comments. 

Remark 2.2.1: In what follows, we suppose that is non trivial. In other 
words, we start with some initial data satisfying 

(2.4) 3(x,e) G M'^ X T; Ul{{),x,e) = UlQ{x,e) ^ . ^ 
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Remark 2.2.2: Fix any Z G N*. The Borel's summation process allows to 
take = +cx) in the Theorem 2.1. It yields BKW solutions (u^,pp which 
solve {£) with infinite accuracy 

dtnl + (u^ • V)u^ + Vp^ = 0(£°°) , div = . A 

Remark 2.2.3: Suppose that the function uq G is a global solution of 

Eulcr equations. Suppose also that the phase 999 £ VV^ is subjected to (2.1) 
on the strip [0, cx)[xIR'^ and that it is a global solution of the eiconal equation 
(ei). Then the Theorem 2.1 can be applied with any T G ■ It means that 
no blow up occurs at the level of the equations yielding the profiles U^, 
and the phases (p^. Yet, non linear effects are present. A 

Remark 2.2.4: The characteristic curves of the field dt + UQ- Vx are obtained 
by solving the differential equation 

dt T{t, x) = Uo{t, T{t, x)) , r(0, x) = x. 

Suppose that the oscillations of the profiles U^q are concentrated in some 
domain D cW^. In other words 

supp^[/*oCD, yke{l,--- ,N}. 
The BKW analysis reveals that for all t e [0, T] we have 

supp^i7*(t,-) C {Tit,x); xeD}, Vfce {I,--- , A^} . 
The phenomena under study have a finite speed of propagation. A 

Remark 2.2.5: The influence of dissipation terms will be taken into account 
in the subsection 4.1. The viscosity we will incorporate is anisotropic. It 
is small enough in the direction V^l in order to be compatible with the 
propagation of oscillations. A 

2.3 End of the proof of Theorem 2.1. 

The Theorem 2.1 is a consequence of the Proposition 4.1 which will be 
stated and demonstrated in the subsection 4.2. Below, we just explain how 
to deduce the Theorem 2.1 from the Proposition 4.1 applied with v = Q. 

• Dictionary between the profiles. Select arbitrary initial data for 

no(O,a;)[7*(O,x,0) G (C/fc)(0, x) G , 1 < A; < iV , 

and arbitrary initial data for 

<^fc(0,x) G l<k<l-l. 
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On the contrary, impose 

(2.5) ipk{o,-) = o, yke{l,--- ,N}. 

The Proposition 4.1 provides with finite sequences 

{Uk}l<k<N , {Pk}l<k<N , Wk}l<k<N , 
and source terms 

such that the associated oscillations 

u^(t,x) := uo(t,x) + Y.k=i Uk{t,x,e-^ ^pl{t,x)) , 
t>l{t,x) := Po(t,a;) + X;f=i ^~ Pk{t,x,e-^ (pl{t,x)) , 
i^%t,x) := f\'{t,x,e-'ip'^{t,x)), 
g^(t,x) := gf{t,x,e-^(pl{t,x)) , 

are subjected to 

5iU^+(ii^-V)u^+Vp^ = \^ = Oie"^-') , div Ci^ = = Oie'^-') ■ 

The oscillations and p^ involve the complete phase f^{t,x) which is the 
sum of the geometrical phase (pg{t,x) plus some adjusting phase eip'^{t,x). 
More precisely 

ip^^{t,x) := ipl{t,x)+e ipl{t,x), Vl{t,x) := YJk=l 'Pk{t,x). 

The functions and can also be written in terms of the phase ipg. 
Indeed, there is a unique decomposition 

iV+l N+1 

= Ub + rUb' = uj' + 0{£~) , Pb = Pb + rPb = Pb + 0(£~) , 
involving the representations 

(2.6) ul{t,x) = ul{t,x,e-^ ipl{t,x)) , p^(i, x) = (t, x, e" a;)) 
where the profiles u^{t, x, 0) and p^(i, x, 0) have the form 

ul{t,x,e) = uo(t,a;) + ^f^^ ei Uk{t,x,e), 

Pl{t,x,e) = Poit,x) + Zk=i Pk{t,x,e). 

The transition from to is achieved through the phase shift 
C/fe(t,x,£-Vb) = ^ik(*>a;,£~Vg + + EfeL/+l Vk)- 
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Use the Taylor formula in order to absorb the small term in the right. It 
furnishes the following explicit link between the {Uk,Pk) and the {UkiPk) 

.27^ Uk{t,x,e-ipiit,x)) := Uk{t,x,9) + g'^iUi,--- ,Uk-i){t,x,e), 

^ ■ ^ Pk{t,x,e-ipi{t,x)) := Pk{t,x,e) + g^{Pu--- ,Pk-i){t,x,e). 

The application Q'^ can be put in the form 

GHUi, • • • , Uk-i) := Ep=! d'^Q^iUi, Uk-p) , k e {!,■■■ ,N}. 
The terms Qp are given by 

Gp{Ul, ■ ■ ■ , ilk-p) ■= ji Y^aejk V^l+l+ai X • • • X (fl+i+ap Ua^+i , 

where the sum is taken over the set 

:= {a= (ai, • ■■ ,ap, a^+i) e N*'+^ ; 

0<aj<N-l-l, e {!,■■■ ,p}, 

1 < "p+i <k — p, oiH \-ap + ap+i = k - p} . 

The relation (2.7) and the definition of Q'' imply that 

Uk{t,x) = {Uk){t,x), "ik e {!,■■■ ,N}, VtG[0,T]. 

Therefore, prescribing the initial data for the Uk or the (Uk) amounts to the 
same thing. The condition (2.5) yields 

g^{Ui,--- ,Uk-pKo,x,e) = o, yke{i,--- ,n}. 

Since ipi{0, •) = 0, we have 

Uo{o,x)u;io,x,0) = Uo{o,x)u;{o,x,e), yke{i,--- ,n}. 

It is clearly equivalent to specify the initial data for the Hq or the Hq . 

• The divergence free relation in the variables {t,x). Consider the 
application 

div : H°° — ^ Im (div) C { g G ; g(0) = } . 
We can select some special right inverse. 

Lemma 2.1. There is a linear operator ridiv : Im(div) — > H°° with 

(2.8) div o ridiv g = g , V g G Im (div) . 

For all L > and for all m G N, there is a constant > such that 

(2.9) II ridiv g Wnm < Cm || g ||^„+i+d+, , Vg e Im (div) . 
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Proof of the Lemma 2.1 . Introduce a cut-ofF function if) G C°°(M'^) such that 
(C; m^^] c 5(0,2], (C; V(0 = 1} ^ 5(0,1]- 

For 5 G Im (div) , take the exphcit formula 

ridiv(g) := (/J V5(Vg)(re) + l^l'^ (1 g(e) x ^ • 

Since g(0) = 0, the relation (2.8) is satisfied. For s > |, the injection 
H^W^) ^ L~(M'^) is continuous. It leads to (2.9). 

• The Leray projector in the variables (t, x). Note n(^) the orthogonal 
projector from onto the plane 

^-L := { u G R'* ; u • ^ = } . 

Introduce the closed subspace 

F := (uG L2; divu = 0} C . 

Call P the orthogonal projector from onto F. It corresponds to the Fourier 
multiplier 

Pu = Yl{D,)n := (27r)-i e'-< n(0fi(6 ■ 

The application P is the Leray projector onto the space of divergence free 
vector fields. It is a self-adjoint operator such that 

ker div = ImP, ImV = (ker (div)) = ker P. 
Consider the Cauchy problem 

dtU + Vp = f, divu = 0, u(0,-)=h 
with data f G L|, and h G L^. It leads to the equivalent conditions 

dtU = P{, u(0,-) = Ph, Vp = (Id-P)f. 

Now we come back to the proof of Theorem 2.1. It remains to absorb the 
term G Im(div). To this end, take i- = jj- Define and p^ as in (2.2) 
with the Uk and P^ of (2.7). Introduce 

cuf := ruf - ridivg^ = 0{ef-'^-i) , cp^ := rpf = 0{s^) . 

After substitution in (S), we lose again a power of e. We find 

f^^ = f^^ - (ridivg^ • V) - (u^ • V) ridivg^ 

— Stridivg^ + (ridivg^ • V) ridivg^ = C){s~i ) . 

The Theorem 2.1 looks like classical statements in one phase non linear 

geometric optics except that the phase c/?^ does depend on e. In the next 
chapter, we examine the part of the (pk which make up (pg and (p^. 
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3 The cascade of phases. 

Turbulence and intermittency are topics which represent extremely different 
points of view. Two approaches compete: 

a) The deterministic approach which study the time evolution of flows arising 
in fluid mechanics [1]- [3]- [11]- [12]- [24]. 

b) The statistical approach in which the velocity of the fluid is a random 
variable [14]- [22]. 

The Theorem 2.1 is mainly connected with a). It brings various informations 
related to the propagation of quasi-singularitics. These aspects are detailed 
at first. Then we briefly explain b) and we draw (in the setting of the 
Theorem 2.1) a phenomenological comparison between a) and b). 

3.1 Microstructures. 

The result 2.1 is concerned with the convection of microstructures. It is 
linked with the multiple scale approach of [24] and [3]. In [24] the authors 
look for BKW solutions in the form 

ul{t,x) = uo{t,x) + U*{t,x,e-H,e-^0oit,x)) + 0{e) . 

In the more recent paper [3] , the selected expansion is 

ul{t,x) = uo{t,x) + £5 Ui{t,x,e~^ t,£-^ 0o{t,x)) + 0{e^. 

Both articles [3] and [24] use homogenization techniques. They perform 

computations involving expressions as or u^. Simplifications (supported 
by engineering experiments) are made in order to get effective equations for 
the evolution of (uq, Uq ) or (uq, Ui). 

Consider the simple case of one phase expansions (that is when ipo = ipo is 
a scalar valued function). Reasons why a complete mathematical analysis 
based on or is not available can be drawn from the Theorem 2.1. For 
instance, look at u^. When I = 3, the oscillation involves the same scales 
as u^g since 

e-^ip'g{t,x) = £-1 (po{t,x)-ire~i (pi{t,x) + e'Ti (p2{t,x). 

Now the analogy stops here since in general ipi{t,x) ^ t and ip2{t,x) ^ 0. 
These are geometrical obstructions which prevent to describe the propaga- 
tion by way of u^. The asymptotic expansion is not suitable. 

Analogous arguments concerning will be presented in the paragraph 3.5. 
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3.2 The geometrical phase. 

Let us examine more carefully how the expression ipg is built. Because of 
the condition (2.1), for e small enough, it is still not stationary 

(3.1) 3eo>0; V(pl{t,x)j^O, V (e, x) G ]0, eo] x [0, T] x M'^ . 

In fact, the function (pg comes from the approximate eiconal equation 

dt^l + H ■ V)^l = Oie) 
which is equivalent to 

at(^jfc + uo • V(^fe + E)=o Uk-j-^Vj = 0, VfcG {I,--- ,1-1}. 

The family {u^(i, 2:)}ee]o.i] e— stratified regularity [19] with respect 

to the phase </?|. This is a geometrical information. 



3.3 Closure problems. 

We have explained why appealing only to (po is not sufficient. It turns out 
that BKW computations relying only on the geometrical phase (p| come also 
to nothing. This is a subtle aspect when proving the Theorem 2.1. We lay 
now stress on it. 

For all iV G N*, the application Q defined below is one to one 
g : (W^)^ — > (W^)^ 



\Un J 



{t,x,9) 



Ui 

U2 + GHU1] 



{t,x,9 + ipi{t,x)) . 



Once the Uj or the Uj are known, it is entirely equivalent to use or u^. 
Before the Uj or the Uj have been identified, in particular when performing 
the BKW calculus, it is deeply different to employ or u^. Indeed, there 

is a unique choice of the (pk with I < k < N, which imposes a specific 
hierarchy between the profiles Uk, which makes possible the triangulation of 
the equations obtained by the formal computations. 

In the subsection 2.3, we will perform the BKW analysis with the profiles 
Uk- It yields a sequence of equations 



(3.2) 



X^'iU^,--- ,Uk+i) = 0, l<k<N. 



14 



As usual in non linear geometric optics, this can be rewritten in order to 
find a sequence of well-posed equations 

(3.3) X^{Uk) = J'iUu--- ,Uk-i), l<k<N, 

where the Uk are made of pieces of the Uj. Of course, the equation (3.3) can 
be interpreted in terms of the Uj and then in terms of the Uj. In this second 
step, it requires to implement the phase shift ipi and the transformations 
with 1 < j < k — 1 and 1 < p < j- Now, the BKW analysis reveals that 
(pi or the various coefficients (pi which appear in the definition of such Qp 
do not depend only on (Ui, • • • , Uk) but also on some Ui with i > k. The 
resulting system is therefore under deter mined. Computations involving the 
functions Uj lead to a sequence of equations which are not closed. 

The insertion of the phases with 1 < /c < iV is an elegant way to introduce 
Q. The change of variables Q, though it is a function of (Ui,--- ,Uis[), is 
needed to progress. It allows to get round closure problems. 

3.4 Compensated compactness. 

Dissipation terms can be incorporated in the discussion. In the variables 
{t, x), the addition of some viscosity k is compatible with the propagation of 
oscillations if for instance k = v . There are approximate solutions (u^, pp 
of the Navicr-Stokcs equations. They satisfy (2.2) and 

dtv^l + (u^ • V)u^ + Vp^ = £2 ^^^e ^ ^ div = , 

with ^^ = 0(£°°)- When > 0, Leray's theorem provides with global weak 
solutions (u^,p^)(t,z) of the following Cauchy problem 

f Stu'^ + (u= • V)u= + Vp^ = v Aa^u*^ , div u= = , 
\ u-(0,.) = u^(0,-). 

Suppose now that Uq = 0. Then, we have also the uniform controls 
(3 4) ^'^^^W ^"^"^ ll-^^l' ^^]0,1]} < C < GO, 

sup { I. £2 II ||2^^^^^^ ; £ G ]0, 1] } < C < GO . 

Arguments issued from the theory of compensated compactness [16] can be 
employed to study the sequence {e^"*" u'^jg. In the spirit of [11] or [12], we 
can try to exploit the informations contained in (3.4) and the equation on 
in order to describe the asymptotic behaviour when e goes to zero of the 
functions £~^ u*^. However this approach seems to be not applicable here. 
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Indeed, obvious instabilities occur. The related mechanisms, which induce 
the non linear instability of Euler equations, arc detailed in the paragraph 
5.1. Below, we just give an intuitive idea of what can happen. Use the 
representation involving the phase ip^. The determination of the inter- 
mediate term tpi requires to identify (Ui) and Ui_i- This is a consequence 
of the equations (4.18) and (4.21). 

In view of the formula (2.7), when tpi is modified by an amount of Sipi, the 
quantity Ui{t,x,9) undergoes a perturbation of the same order 5(pi. When 
dealing with quasi-singularities, some quantities with e in factor (like (f//)) 
or with e^~T in factor (like Ui_i) can control informations of size e~ . This 
fact is expressed by the following rules of transformation 

(3 5) ^^'^ I ^^'^ + ^^^'^ =^ / + O(e^) HUi) , 

UU I UU + SUU =^ / + 0{e7) 5U( . 

Now reverse the preceding reasoning. To describe features in the principal 
oscillating term e~ t/j* (t, x, e~-^ (/9^(t, x)) , we must identify tpi which means 
to obtain (JJi) and In other words, we need to know quantities which 

have respectively e and e^~T in factor. When I > 2 such informations are 
clearly not reachable by rough controls as (3.4). 

This discussion indicates that the study of turbulent regimes requires to 
combine at least geometrical aspects, multiphase analysis and high order 
expansions. The tools of non linear geometric optics seem to be appropriate. 
Some attempts in this direction have already been made. 



3.5 Non linear geometric optics. 

We make in this paragraph 3.5 several comments about non linear geometric 
optics. They concern both old [18]-[19]-[26] and recent [6]-[7]-[8] results 
which all are devoted to one phase expansions of the type 

(3.6) uf (t, x) := uo(t, x) + ^^=1 {t, x, e"! ^o(t, x)) . 

When / = 1, one is faced with weakly non linear geometric optics. The 
asymptotic behavior and the stability of are well understood. In fact a 
complete theory has been achieved in the general framework of multidimen- 
sional systems of conservation laws (see [18]- [19] and the related references). 
Because of the formation of shocks, the life span of exact solutions close to 
uf does not go beyond T ~ 1. 
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When 1 = 2, expressions as are called strong oscillations. The construc- 
tion of such BKW solutions can be undertaken only if the system of conser- 
vation laws has a special structure. Transparency conditions are needed to 
progress. They can be deduced from the presence of a linearly degenerate 
field [7]. In the hyperbolic situation the family {u^}ee]o,i] is unstable [7] on 
the interval [0,T]. It becomes stable on condition that a small viscosity is 
incorporated [6]. Applications can be given to describe large-scale motions 
in the atmosphere [6]. 

Compressible Euler equations are the prototype of a non linear hyperbolic 
system having a linearly degenerate field. After a finite time, singularities 
appear. These correspond to the generation of shocks by compression [27]. 
The situation is different in the incompressible setting. There is no genuine 
shock and the production of singularities poses a much more subtle problem 
[2]- [9] which up to now remains basically open. 

Incompressible fluid equations lie at an extreme end in the sense that they 
are the most degenerate (or the most linear) equations which have just been 
mentioned. Following the approach of [20] related to transparency, repeating 
the reasoning which goes from [18]- [19] to [6]- [7], one expects to go further 
than / = 2 when dealing with {£). Now, this is precisely what says the 
Theorem 2.1 since it allows to reach any Z G N* ! 

To tackle the limit case Z = oo, one is tempted to look at asymptotic expan- 
sions of the form 

(3.7) nl,{t,x) := Ek=o^'' Uk{t,x,e-\o{t,x)) , deU^ ^ . 

The oscillations contained in have a large amplitude. Modulation equa- 
tions for Uq are proposed in [26]. However these transport equations are 
not hyperbolic so that they are ill posed (in the sense of Hadamard) with 
respect to the initial value problem. It confirms that a BKW construction 
based on (3.7) is not relevant^. 

The contribution [26] does not explain why the expansion (3.7) is not the 
good one. We come back below to this point. At first sight the Theorem 
2.1 does not include large amplitude waves since — uq = 0(^^) 0(1)- 
A change of variables leads to recant this impression. Suppose that uq = 
and deU^ # 0. Then define 

u;(/..r) := c^T (.--T /. ,x-) , pt{t,x) := e'T p^{e-Jt,x). 

^The singularities are carried here by the velocity field. The discussion is very different 
when the oscillations are polarized on the entropy [8]. 
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Observe that the structure of and is very different from the one in 
(3.7) since we have 

Ub(*'-^) = T.kLi Uk{e-T t,x,e-^ iflie-T t,x)) +e"T cu^(e"T i, x) , 
Pb(*'^) = Efcli Pk{e-h, X, e-^ iflie'T t,x)) +£-? cpl{£~Jt,x). 

The functions li^ and satisfy 

atU^+(u^-V)u^+Vpf = divu^ = 0, f^%t,x) = e-J {^%e-Tt,x). 

The functions are oscillations of the order 1. They are approximate 

solutions of {£) on the small interval [0,e~ T]. Indeed, for all m G N, the 
family {(^}e is subjected to the uniform majoration 

iV_j_ 2_i_Q_i_ 

SUp,g]o,,„] £-T + T+3+- II fe \\^rn^^^^ ^ < OO. 

If moreover TV = +00 and 

MO,-) = ■■■ = v^i-iio,-) = 0, 

^""•^^ C/fe+i(0,-) = 0, VA;gN\(ZN), 

the trace u^(0, •) has the form 

Ub(0'-^) = Er=o f/i+/fe(0,x,e-Vo(0,x)) , a^f/r ^ 0. 

At the time t = 0, we recover (3.7). Now the construction underlying the 
Theorem 2.1 reveals that in general 

(3.9) ipk{t,-)^0, VtG]0,r], VfeG{2,--- ,/-l}. 

The functions tpj with j G {2, • • • ,1 — 1} are not present when t = 0. But 

1 k 

the description of u^(t, •) on the interval [0, e ~~ T] with A; G {2, • • • , Z — 1} 
requires the introduction of the phase shifts ipj for j G {2, • • • ,k}. More 
generally, the description of u^(t,-) on the whole interval [0, T] needs the 
introduction of an infinite cascade of phases {<fj}jeN*- 

Such a phenomenon does not occur when constructing large amplitude oscil- 
lations for systems of conservation laws in one space dimension [10]- [13]. It 
is specific to the multidimensional framework. It explains why the classical 
approach of [26] fails. 

It seems that the creation of the ipj is due to mechanisms which have not 
already been studied. It is not linked with resonances. It is related neither 
to dispersive nor to diffractive effects. 
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Remark 3.5.1 (about ipi): The term ipi does not appear if (pi{0,-) = and 
Ui{0, ■) = 0. When these two conditions are not verified, the phase shift ipi 
can be absorbed by the technical trick exposed in [7]. Just replace uo(0, •) 
by Uo(0, ■) + 6 Ui{0, •). Perform the BKW calculus with a fixed 6 > 0. Then 
choose S = e. A 

Remark 3.5.2 (about ip2)-' In general, we have <f2^0 even if 

MO, •) = ^2(0, •) = , Ui{o, •) = [72(0, •) = . 

Indeed the time evolution of U2 is governed by (4.22). It involves the source 
term div (?7j* Ul) which is able to awake U2. This infiuence can then be 
transmitted to (/72 through the transport equation 

(3.10) dtip2 + (uo • V)ip2 + {Ui ■ V)</?i + {U2 ■ V)^o = . 

Likewise, the other terms (ps, • • • , <^/-i are in general non trivial even if 

(^1(0, •) = ••• = fi-iiO, ■) = 0, [7i(0, .) = ... = Ui^O, .) = . 
There is no more trick which allows to get rid of (^2 j • • • , <fi-i- A 

Remark 3.5.3 (why turbulent flows ?): The introduction of the phase shifts 
(pk with 2 < k < I — 1 cannot be avoided. Therefore the difficulties that we 
deal with appear from I = 3. When I > 3, the characteristic rate e of eddy 
dissipation is bigger than one [3] . This is the reason why such situations arc 
refered to turbulent regimes. A 

Remark 3.5.4 (about shear layers): We have said in the introduction that 
the expression u| given by formula (1.1) is of a very special form. Let us 
explain why. Change the variable t into eT t and into iif := u^. The 
main phase LpQ(t,x) = X2 remains the same. Now we are faced with 

mt,x) := *{eT g{x2,e-^ X2),0,eT h{xi-eT g{x2,e-^ X2)t,X2,e-^ X2)) . 

It is still a solution of Euler equations. Now it falls in the framework of the 
Theorem 2.1. The constraints on tJ2 = ^([^2 ' ^2 ' ^2 ) reduce to 

C/i = C7| = 0, dttJl + {gdih) = 0. 

The contribution U2 is non trivial but it is polarized so that tJ2 ■ Vi^o = 0. 
Therefore it does not produce the phase shift (p2. The same phenomenon 
occurs concerning ip^, ■ ■ ■ , ipi-i. These terms are not present. It turns out 
that the expansion involves only the phase ipQ{t,x) = X2. A 
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The choice for the ampUtude of the oscillations is very important. It is 
strongly related to the scale of time T under consideration. The idea is to 

increase the time of propagation T to reach the regime where non linear 
effects appear. Starting with some large amplitude high frequency waves 

u^(0,x) = f/o(0,,x,e-Vo(0,.x)) + Qis) , deUS{0,-)^0, 

the preceding discussion can be summarized by the following diagram: 



r~ 1 

infinite cascade 

of phases 
<P0 - {<Pl) 

VO - (Vl) - 

r ~ £^ 

T c:±e 

r = o 

phases 



turbulent 
flows 



turbulent 
flows 



strong 
oscillations 
[6] - [7] 

weakly 
non linear 
geometric 

optics 
[18] - [19] 

regimes 



incompressible 
fluid equations 



incompressible 
fluid equations 



systems of conservation 
laws with a linearly 
degenerate field 



systems 
of 

conservation 
laws 



equations 



This picture allows to understand the position of the actual paper in com- 
parison with previous results. 

3.6 The statistical approach. 

It deals mainly with quantitative informations obtained at the level of ex- 
pressions, say u(a;), which in general do not depend on the time t. The 
introduction of u can be achieved by looking at stationary statistical solu- 
tions [14] of the Navier-Stokes equations that is 

u{x) = limy >oo 7 /o^ u{t,x) dt 
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or in conjunction with the ensemble average operator ([22]-V-6) marked by 

the brackets < • >. We will follow this second option. The description below 
is extracted from the book of M. Lesieur [22] (chapters V and VI). We work 
with d = 3. Interesting quantities are the mean kinetic energy 

i <U(X)2 > ~ ^3 |u(.t)|2 dx, 

the enstrophy (that is the space integral of the square norm of the vorticity) 

i < u;(a;)^ > ~ f^^ \uj{x)\'^ dx , uj{x) ■.= V Au{x) 

and the rate of dissipation e ~ k < lo{x)'^ >. In the setting of isotropic 
turbulence, these quantities can be expressed in terms of a scalar function 
k I — > E{k). The real number E{k) represents the density of kinetic energy 
at wave number k (or the kinetic energy in Fourier space integrated on a 
sphere of radius k) . The relations are the following 

[22]-V-10-4 i < u(x)2 > = /+°° E{k) dk . 

[22]-V-10-15 i < a;(x)2 > = /+°° P E{k) dk . 

[22]-VI-3-15 e = 2k /+°° k"^ E{k) dk . 

Kolmogorov's theory assumes that 

[22]-VI-4-l 3 c> ; E{k) = c e^/^ A;-5/3 ^ V A; G [iti , . 
This law is valid up to the frequency kd with 
[22]-VI-4-2 kd ~ (e/K3)V4. 

The small quantity s := kj^ is the Kolmogorov dissipative scale. The rela- 
tions [22]-VI-3-15 and [22]-VI-4-2 imply that the rate of injection of kinetic 
energy e is linked to the number I according to e ~ £ ' . We recover here 
that e ~ 1 when / = 3 (see [3]). 

A starting point for the conventional theory of turbulence is the notion that, 
on average, kinetic energy is transfered from low wave numbers modes to 
high wave numbers modes. A recent paper [14] put forward the following 
idea: in the spectral region below that of injection of energy, an inverse 
(from high to low modes) transfer of energy takes place. At any rate, it is a 
central question to determine how the kinetic energy is distributed. 

3.7 Phenomenological comparison. 

The statistical approach is concerned with the spectral properties of solu- 
tions. Below, we draw a parallel with the propagation of quasi-singularities 
as it is described in the Theorem 2.1. 
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Let us examine how the square J^{u^){t,^)^ of the Fourier transform of 
u^{t,x) is distributed. To this end, consider the apphcation 

E{t, ■) : R+ — ^ M+ 

The initial data u^(0, •) has a spectral gap. In another words, the graph of the 
function k i — > E{0, k) appears concentrated around the two characteristic 
wave numbers A; ~ 1 and k = k^- In view of (3.9), this situation does 

not persist. At the time t = , the concentration is around I characteristic 
wave numbers which are intermediate between the two preceding ones. This 
corresponds to a discrete cascade of energy. 

Suppose now (3.8) and consider li^. The life span of u^(t, •) is ei T. There 
are various manners to get a family {uf (t, OleeJO,!] which is defined on some 
interval [0, f] with T > independent on e. In particular, we can 

a) Select any T > when T = +oo. However nothing guarantees that the 
functions li^ are still approximate solutions on the interval [0,T]. Indeed, 
since t is replaced by £~(^/') t, the size of the error terms depends on the 
increase of with respect to t. At this level, we are faced with secular 
growth problems [21]. 

b) Use a convergence process^ which needs the introduction of an infinite 
cascade of phase shifts. The intuition^ is that the graph of E becomes 
continuous (no more gap). This corresponds to the impression of an infinite 
cascade of energy. This remark is consistent with engineering experiments 
and the observations reported in the statistical approach. 

The turbulent phenomena which we study are very complex in their 
realization. When t > 0, the description of u^(i, •) involves an infinite set 
of phases so that computations and representations are hard to implement. 
It gives the impression of a chaos. Nevertheless, our analysis reveals that 
these phenomena contain no mystery in their generation. On the contrary 
quantitative and qualitative features can be predicted in the framework of 
non linear geometric optics. 



^When performing the formal analysis, arbitrary values can be given to the parameters 
£ € ]0, 1] and / G N* . For instance e can be fixed whereas I goes to oo. Or / = — (In e)/ (In 2) 
so that eT T = iT > 0. 

''Even at a formal level, difficulties occur in order to justify the different convergences. 
Rigorous results in this direction seem to be a difficult task. 
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4 Euler equations in the variables {t,x,9). 



As explained in the previous chapter, the demonstration of the Theorem 2.1 
is achieved with the representation 

(4.1) mt,x) =ul{t,X,£-'^ipl{t,x)) , p^^{t,x) =p^^{t,X,£-^(p^^{t,x)) . 

Recall that the complete phase (p^{t,x) is 

(4.2) (pl{t, x) = (plit, x)+e iplit, x) = (po{t, x) + Ylk^i £^ <^fe(i, x) 
and that the profiles «^(t, x, 6) and p^(t, x, 9) have the form 



ul{t,x,9) = uo(t,x) + XlfcLi Uk{t,x 



(4.3) 



Pl{t,x,9) = Po(i,x) + Ef=i £^ Pk{t,x,9). 
4.1 Preliminciries. 

• Anisotropic viscosity. Mark the abbreviated notations 

X^it,x) :=V(^f(t,x) = Ef=o^^ Xk{t,x), Xk{t,x) :=V^fe(i,x), 
Xl,{t,x) ■.= \X^it,x)\-^ X^{t,x). 
Complete the unit vector X^i{t, x) into some orthonormal basis of M"* 



Xl^{t,x)-Xl^it,x) = 6, 



V(i,i) G {!,••• ,d}' 



The 



so that all the vector fields X^- are smooth functions on [0, T] x I 
corresponding differential operators are denoted 

^Ud):=^Ut,x)-V, ie{l,---,d}. 
Their adjoints are 

^iid)* :=X^,(t,a;)-V + div(X^.)(t,x), i e {!,■■■ ,d}. 
Select q G C^{[0,T] x R"^; Sf) be such that 

3c>0; q{t,x)>c, V (t, x) G [0, T] x M'^ . 
Let {m,n) G N^. Consider the elliptic operator E^^{d) defined according to 

Eern^QJ := {sf ^,{8}* , ST ^^,(5)* , ■ ■ ■ , sf X^,(9)* ) 



/ qii(t,x) qi2it,x) 
(\2i{t,x) q22{t,x) 

\ c\di{t,x) qd2{t,x) 



^ld{t,x) \ 

<\2dit,x) 



e" XUd) 



c\dd{t,x) J V s'^XfM J 
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The introduction of the operator E^^{d) in the right of {£) is compatible 

with the propagation of oscillations only if m > / and n > 0. We retain the 
limit case / = m and n = 0. The other situations are easier to deal with, at 
least when performing formal computations. 

• Interpretation in {t,x,9). To deal with the variables {t,x,6), define 
:= £ dj + djifl xde, j e {0, - ■ ■ ,d} , 

Qtadl := \di,e, ■ ■ ■ ,T)d,e) = sV + X^xde, 

t)it)^ := (sroDp* = £ div + ■ de . 

The derivatives become 

eXl^{^e) -e^^id) + \X'^{t,x)\xde, 

£X^.(c),) ■.= eXl^{d), Vj€{2,... ,d]. 

The action of E^^{d) expressed in the variables (t, x, 6) gives rise to some neg- 
ative differential operator of the order two, noted E^^{ps). The coefficients of 
the derivatives in E^q{T)^) are of size one, except in front of Xl^{d). To avoid 
technicalities and to simplify the notations, we substitute the Laplacian v A 
for Ell^ili,). 

When z/ = 0, we recover Euler equations. When z/ > 0, the action z^A can 
be viewed as the 'trace' in {t,x,6) of the anisotropic viscosity £'^q(9). Now, 
consider the Cauchy problem 

( u^(0,-) = hli-), 
with given data 

Suppose that = and select some smooth solution {u^,pl) of (4.4). The 
expressions and given by the formula (4.1) are subjected to 

dtui + (li^ • V)u^ + Vp^ = f^^ , div nl = gl, 

u^(o,•) = Hi-), 



(4.4) 



(4.5) 



where the functions f^{t,x), g^(t,x) and h^{t,x) are obtained by replacing 
the variable 9 by ^pl(t, x) in the expressions f^i^^ 2;, 9), ffb ^) 
h^{t, X, 9). In other words, any solution of (4.4) with v = Q yields a solution 
of (4.5). Prom now on, we proceed directly with the relaxed system (4.4). 
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4.2 The BKW analysis. 

Select a smooth solution uo(t, x) G >V|P of 

dtuo + (uo ■ V) uo + Vpo = Aa;Uo , div uo = . 

Choose a phase (po{t, x) G Ci([0, T] x M'^) with Vv^ol*, x) G C^°°([0, T] x R"^). 
Suppose moreover that it satisfies the eiconal equation (ei) and the condition 
(2.1). The main step in the construction of approximate solutions is the 
following intermediate result. 

Proposition 4.1. Select any b = {I, N) G such that < I < N. Consider 
the following initial data 

U*f,{x, 9) = no(0, x) C/*o(x, e)eH°°, l<k<N, 

{Uko){x)eH^, l<k<N, 
(Pko{x)eH'^, l<k<N. 

There are finite sequences {Uk}i<k<N o-'f'd {Pk\i<k<N with 

Uk{t,x,9) eW^ , Pk{t,x,e) eW^ , l<k<N, 
and a finite sequence {<fk}i<k<N with 

(Pkit,x)eW^, l<k<N, 
which are such that 

no(0, x) U*{0, X, 9) = no(0, x) UIq{x,9), l<k<N, 

{Uk){0,x) = {Uko){x), l<k<N, 

ipkiO,x) = (pko{x) , l<k<N. 

Define (p^ as in (4-2)- All the preceding expressions are adjusted so that the 

functions and associated with the expansions in (4-3) are approximate 
solutions on the interval [0,r]. More precisely, they satisfy (4-4) with 

(4.6) hl{x,9) = uo(0,x) + Ef=i^^ Uko{x,9) 
and we have 

(4.7) {a = 0(£^), m. = 0(£^)- 

• Proof of the Proposition 4.1. For convenience, we will drop in this 
paragraph the tilde '~' on the profiles u^, p^, Uk and P^. This modifica- 
tion concerns only this demonstration. We hope that it will not induce 
confusions: we still work here with the complete phase (p^. 
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Because of (3.1) we can define the application Il^{t,x) which is the orthogo- 
nal projector on the hyperplane V(p^{t, x)-^ C W^. We adopt the convention 

The access to 11^ needs only the knowledge of the Xj for j < k. Introduce 



-.k-l 



By construction 

< = «b l^b r' + < ' Uk = Vk |Xo|-2 Xo+Wk + h 
where h depends only on the Xj for j < k and on the Uj for j < k — 1. 

The conditions prescribed in the Proposition 4.1 on the initial data U^q allow 
to fix the functions V(/?o(0,x) • U^{0,x,6) as we want. Since 

V* = V^o • u* + Xj . u*_. , 

the same is true (by induction) for the components V^{0,x,6). To begin 
with, we impose the polarization conditions 

(4.8) P*^v*^o, yke{l,---,l} 

and we adjust a priori the geometrical phase ipg so that 

(4.9) dtiPk + Vk = o, yke{i,--- ,1-1} 

which implies that 

dtifl + (n^ • V)ipl = Ek=i Vk = 0(e) • 

It amounts to the same thing to look at the equations in (4.4) or at the 
following singular system (we drop here the indices e and b at the level of 

■f^b' pI ' ipl) 



(4.10) 



dtu + {u ■ V) u + Vp + [dfif + v) dgu 

+ £~'^ dep V(/? = V /S.U, 

div u + dgv = . 



The functions v is subjected to 

(A-i-i) dfV + {u ■ V) V + X ■ Vp + {dt^ + v) dev 

^ ' +e-^ dep \\Xf-{dtX + {u-V)X)-u = vX-I!^u. 
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The functions w satisfies 



(4.12) 



dtw + {u ■ V) w + U Vp + e ^ {dt(p + v) dew 



Substitute the expressions and given by (4.3) into (4.10). Then arrange 
the terms according to the different powers of e which are in factor. The 
contributions coming from the orders £~~^, • • • , £~~ and e*^ are eUminated 
through (4.8), (4.9) and the constraints imposed on (uo,po)- 



r OtU, + Ei=n {Uk • V) Uj.k + VPj + Ei=o dePi+k y^j-k 
(4.13) \ + Y:i-Jo {dm+k + Vi+k) deUj^k = iyAUj, 



Proceed in a similar manner with (4.11). Just arrange the terms which come 
from X ■ Au and which do not involve Uj in a source term 



dtVj + n=o (Uk • V) Vj.k + E£o (dm+k + Vi+k) deVj.k 

+ Ei=o ■ vPo-k - Ei=o dtXt ■ u,-k 

(4-14) - Ei=o (Eto {^^-i ■ V) X{) . C/,_fe 

+ EU (Ezio ■ Xk-i) dePj+i-k + |Xop dePj+i 

= uXo-AUj+ n'yit, X, ^, f/i, Xi, • • • , f/,_i, X,) . 



dtWj + Ei.=o iUk ■ V) + ECo (^tV^'+fc + 9eWj.k 

- ELo (Ef=o iu,-i • V) no 

= uUoAUj + n\^{t, X, 0, C/i, Xi, • • • , Xj) . 



Now, look at the terms in front of ei with j G N*. It remains 




The same operation with (4.12) yields 



Then extract the mean value of (4.13) 



r dtUj + (uo • V) tjj + {Uj • V) uo + VP,- 



(4.16) 



+ ECl((^^-v)t/,-fc) 

+ Y.'kMyi\kdeUj-k) =i^A,Uj, 
, div Uj = 0. 



Observe also that 



(4.17) 



-div dn^U* 



3 ' 
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Consider the inductive reasoning based on 
Hypothesis (Hj) : 

i) The expressions Ui, ■ ■ ■ , Uj and Pi, ■ ■ ■ , Pj are known. 

ii) The phases (pi, ■ ■ ■ , ipj are identified. The same is true for the vectors 
Xi, ■ ■ ■ , Xj and the projectors Hi, ■ ■ ■ , Uj. Moreover, the following relations 
are satisfied 

(4.18) dtiPj+k + Vj+k = 0, VfcG {!,••• ,^-1}. 

Hi) The correctors V^+i, • • • ; o,i^d -P^+i, • • • , Pj^i are identified and 

(4.19) = - div d^^u;^k_i , yke {!,■■■ , I}. 

- Verification of (Hi). The mean value C7i is obtained by solving 

dtUi + (uo ■V)Ui + (Ui ■ V) uo + VPi = iyA^Ui, div t/i = . 

Using (4.9), it allows to determine ipi. Now look at the oscillating part of 
(4.15) with the indice j = 1. The constraint on Wf = writes 

dtW* + (uo • V) w* + (dm + Vi) deW* = mw* + i^Uo a 

where M is the linear application 

MU := (dtUo) U + ((uo ■ V)no) U - Uo {U ■ V)uo . 
We impose (4.18) for j = 1. In view of (4.9), it reduces to 
dtipi + Vi = 0. 

The link between and Vi is removed. It remains the linear equation 

dtW^ + (uo • V) = MWl + uUo A . 

When u = 0, the profile W* is obtained by integrating along the character- 
istic curves r{-,x). This justifies the remark 2.2.4 for A; = 1. Observe also 
that the polarization condition = Uq is conserved since the equation 
given for Wf is equivalent to 



(4.20) 



no [dtWf + (uo • V) + {Wf ■ V)uo] = uIloAWf, 



Introduce the linear form 

iU := |Xo|-2 [dtXo ■ U + ((uo • V)Xo) ■ U - Xo ■ {{U ■ V)uo) ] . 
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The constraint = is equivalent to 

PUi = e dg^W^ + u |Xo|-2 Xo ■ dg^ A , 
We have also 

At this stage, we know who is Ui = Ui + Wi and Pi = 0. Moreover, we 
have the relations (4.18) and (4.19). Thus, the hypothesis (Hi) is verified. 

- The induction . Suppose that the conditions given in (Hj) are satisfied. 
The question is to obtain (iJj+i). Consider first (4.16) with the indice j + 1. 
The relation (4.19) induces simplifications. It remains 



(4.21) 



( dtUj+i + (uo • V) Lfj+i + {Uj+i • V) uo 

+ Ei=i (Uk ■ V) Uj+^-k + Ei=i div {u* ® 

+ VPj+i = u Uj+i , div Uj+i = . 



This system gives access to Uj+i and Pj+i- For j = 1, it yields 
(4.22) 



dtU2 + (uo • V) C/2 + {U2 • V) Uo + VP2 

+ (Ui ■V)Ui + div {UI (g) UI) = V AU2 , div f/2 = . 

Because of (2.4), the source term {U^ Uf) is sure to be non trivial. We 
recover here that in general U2 ^ even if Ui{0, •) = L^2(0, •) = 0. The term 
U2 excites ip2 through (3.10). Generically, we have (^2^0 even if 

c7i(o,-) = U2{o,-) = 0, (^1(0, •) = ^2(0,0 = 0. 

Observe however that exceptions can happen (see the remark 3.5.4). The 
information (4.18) for A; = 1 means that 

dtfj+i + (uo • V) ipj+i + Xo ■ Uj+i + X;/=i Xi ■ Uj+i^i = . 

Deduce ipj+i from this equation, and therefore Xj^i and Ilj+i. Complete 
with the triangulation condition 

(4.23) dt^j+i + Vj+i = 0. 

Then extract the oscillating part of (4.15) written with j + 1. Use (Hj) and 

(4.23) in order to simplify the resulting equation. It yields 

(4.24) dtW*+^ + (uo • V) W^+i = M t^jVi + iyUoA W*+^ + f 

where / is known. We get Wj_^_i by solving (4.24). Therefore we have t^^+i 
and we can deduce V^+j+i = — div 9^^f7*_|_^. 
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Now look at the constraint (4.14) for the indice j + Extract the oscillating 
part. It allows to recover Pj_^i^i. Thus we have {Hj^i). 

Apply the induction up to j = — /. It yields Ui, ■ ■ ■ , Un- Construct 
oscillations and by way of (4.3). It furnishes source terms and 
through (4.4). By construction, we recover (4.7). □ 

4.3 Divergence free approximate solutions in {t,x,9). 

In this subsection 4.3, we impose on ipo a constraint which is more restrictive 
than (2.1). We suppose that we can find a direction ( eW^X {0} such that 

(4.25) 3c>0; V^o{t,x) ■ ( > c, V (t, x) G [0, T] x R'^ . 

Proposition 4.2. The assumptions are as in the Proposition 4-1 ■ The 
profiles Uk, Pk> o-n-d the phases (pk are defined in the same way. Then, there 
are correctors 

culit, X, 9)eW^, {cul), = Oief) 

such that the functions and pf defined according to 

ul{t,x) := uo{t,x) + Ylk=i Uk{t,x,e~'^ ipl{t,x)) +cul{t,x) 
pl{t,x) := Po(i,a:) + X)feLi ^"^ A;(*, a;, e" Vg(i, a;)) 

satisfy the Cauchy problem 



(4.26) 



. ui{o,x,e) = uoio,x) + Y:k=i £^ Uko{x,e) 



and we still have {f^}e = 0(^~'~)- 

We need some material before proving the Proposition 4.2. 

• The divergence free relation in the variables (t,x,9). We can select 
some special right inverse of the application Dw^ : H^* — > H^*. 

Lemma 4.1. There is a linear operator riOiti^ : Im(t)it)p — > H^* with 

(4.27) X)it)l o ritJlD^ 9 = 9, V 5 € Im (310^ . 
For all m eN, there is a constant Cm > such that 

(4.28) II tidiol 9\\Hm< Cm \\ 9 II ^^+1+4 , G Im (Oiep . 
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Proof of the Lemma 4^- Let n G N*. Note 

tj := jT/n , Xj = k/n , 1<J<?^ — 1, /cGZ*^. 
Consider a related partition of unity 

Xim G C'-([0, T] X M'^) , (j, fc) G {1, • • • , n - 1} X 

E;=i' Efcez^ X0-,fc)(t, x) = 1 , V (t, x) G [0, T] X 
{ (i, ^) ; X(i,fc)(i, / } C [t, - f , + f ] X f ] , 

{ {t, x) ; X(i,fc)(i, x) = l} D [t^- - + i] X i] . 

By hypothesis, there is a function v G such that g = t)iD^ Introduce 

It suffices to exhibit ritiiti^ 9(j,k) to show (4.28) with a constant Cm which 
is uniform in {j,k). The problem of finding tidit^ 9ij,k) can be reduced to 
a model situation. This can be achieved by using a change of variables in 
(t,x), based on (4.25). From now on, the time t is viewed as a parameter, 
the space variable is x = (a;i,x) G M X W^'^, and we work with 

g = g* = v = {edi + do)vi + d2V2 H h dava , 

(x; g{x,e)^0} C {x; v{x,e)^0} C B(0, i] . 

Let G C°°(M'^-^;R+) be such that J^^-i '^{x) dx = 1 and 

{x; 4^{x)^0} C 5(0,1], (f; V(*) = 1 } D 5(0, . 

Decompose 5 according to 

9 = {9-9)'<P + 9'^, 9{x) ■■= /iKd-i g{xi,x) dx = (e^i + d0)vi . 

Seek a special solution u having the form 

u = tidit)! g = \a, ridiv [{g - g) i)]) , a G Hf* 

where 'ridiv' is the operator of Lemma 2.1 applied in the dimension d — \. 
It remains to control the scalar function a which satisfies the constraint 

e d\a + dea = h := g ip = {e d\ + de){v\ tp) . 
Take the explicit solution 

a(xi, X, 6) = /i(xi + £ (s — 9), X, s) ds 

= /°oo ^(^1 + r,x,9 + r) dr . 
By construction 

a{x,9 + 1) = a{x,e) , J,^ a{x,9) d9 = , V (x, 6*) G M'^ x T . 
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For |a;i| + \x\ > 2, we find 

(^(^,0) = i [ivii^){xi + e{s-9),x,s)] ds = {hi^)ix,e) = 0. 
It implies that 

{ {x,e); a{x,e) /O} C 5(0; 2]. 
Note f) := dg^h G Obviously 

II II//™ < Cm II h ||i/rn , Vm G N, 

{{x,9); i){x,e)^0} C 5(0; 1], 
and we have the identity 

a{xi,x,e) = t){xi,x,e) - JI^^^^ di\:){xi+r,x,9 + £-^r) dr. 

The term on the right is supported in 5(0,2]. Use Fubini and Cauchy- 
Schwarz inequality to control the integration of dii). It yields (4.28). 

• The Leray projector interpreted in the variables {t,x,6). Introduce 
the closed subspace 

F^ := {u* e ; VW^u* = } C L^* . 

Note the orthogonal projector from L^* onto F^. This is a self-adjoint 
operator such that 

ker Oit)^ = Im^}^ , Imgrot)^ = (ker Oiop"^ = ker . 

Expand the function u* € L^? in Fourier series and decompose the action of 
in view of the Fourier modes 

u*{t,x,e) = ^,^^^ u,(^,x)e'^^ ?Jfn* = Eik6Z. VlkMt,x) e^'^r 
Simple computations indicate that 

The following result explains why the projector is replaced by IIq when 
performing the BKW calculus. 

Lemma 4.2. 

i) The family {^l}e is in We have [de;%] = and 

a) The projector Il^{t, x) is an approximation of^^ in the sense that 
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Proof of the Lemma 4-2. Since is a projector, we are sure that 

II « - II ^ II^T ' ^ ^) ^ ^ ■ 

It shows that {^l}e G Compute 

Observe that 

(£ dj + ik djipl) qj^^Ufc = e-'^-'^'^t U{D^) edj {e''~'^^l u^) 

= %k dj + ik djipl)uk. 

All these informations give access to the first assertion i). Now consider ii). 
The asymptotic expansion formula for pseudo differential operators say that 
for all ujfc in C^(M^) we have 

V(i,x)GM^, ^lim^ {{%,^nk){t,x) - Ii{y^l{t,x))uk{t,x)] = 0. 

Since 11^ = n(V(/?p, it indicates that is close to 11^. Wc have to make 
this information more precise. To this end, proceed to the decomposition 

u* = V* + e Vp* + dep* X , v* = ^lu* . 
We seek a solution {v*,p*) of these constraints such that 

V* = U^u* + e V* , p* = II ||-2 X^ ■ dg'^u* + ep*. 
After substitution, we find the relation 

- V(|| X^ II -2 X^ ■ dg^*) =v* + eVp* + dep* X X^ 
which must be completed by the condition 

- div (n^ u*) = e divv* + X^ ■ dev* . 
It follows that 

V* = -% [v(|| x^ ||-2 x^ . d^\*)] + m - Id) tioio^ (div (n^- u*)) . 

In view of this relation, the point ii) becomes clear. 
Consider the Cauchy problem 

X)o,eU* + e-^ Qtad^p* = r, X)Wlu* = 0, u*{0,-) = h*{-) 
with data /* G Lf? and h* G L^*. Compose on the left with It yields 

The Cauchy problem can be solved in two steps. First extract u* from the 
above equation. Then recover p* from the remaining relations. 
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• Proof of the Proposition 4.2. It remains to absorb the term g^. Use 
the decomposition 

9l = {qD + ~9l* , {9l) e Im (div) , ~gl* G Im (Oiop . 
It suffices to choose 

AT+l 

cuf := - ridiv {gD - tiOiuf gl* = 0(£~) • 

5 Stability of strong oscillations 

The case of turbulent regimes (Z > 3) will not be undertaken here. Prom 
now on, fix Z = 2 and AT » (6 + d). Consider the Cauchy problem 

, . f dtu' + (u^ • V) + Vp= = u Ef^^{d) , div = , 
^^•-"^ \ u-(0,x) = u^(0,x). 

Let be the upper bound of the T > such that (5.1) has a solution 
G Wji. Classical results [4] for fiuid equations imply that > 0. Our 

aim in this chapter 5 is to investigate the singular limit 'e goes to zero\ Such 
an analysis must at least contain the two following parts. 

a) An existence result for a time Tq which is independent on the small 
parameter e G ]0, Eq]. It is required that 

inf {T,; eG]0,l]} > To > 0. 

When v > 0, OT when z/ = and d = 2, we know [4]- [23] that Tg = +oo so 
that To = +00. When i/ = and d> 3, nothing guarantees that Tq > 0. To 
our knowledge, this is an open question. 

b) A convergence result. The exact solution is not sure to remain close 

on the whole interval [0, Tq] to the approximate solution given by the 
Theorem 2.1. Proving estimates on u*^ — is a delicate matter. 

5.1 Various types of instabilities 

• Obvious instabilities. The obvious instabilities are the mechanisms 
of amplifications which can be detected by looking directly at the formal 

expansions u^. They imply the non linear instability of Euler equations. 
Indeed, fix any T > 0, any uq G W^{W^) which is solution of {£), and any 
(5 > 0. Work on the balls 

5o(uo; 5] := { u G L2 ; II u(-) - uo(0, •) ||i2(Rd) < S} . 
Bt{uo;6] := { u G L|, ; || u - Uq ||L2([o,r]xM<*) < S} . 
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Proposition 5.1. For all constant C > 0, there are small data 
(h,h) G {Bo{uo;5]nH°°)\ (f,f) G {Bt{uo;6] nW^f 

so that the Cauchy problems 

atu + (u- V)u + Vp = f , divu = 0, u(0,-) = h(-), 
atu + (u- V)a + Vp = f , divu = 0, u(0, •) = h(-)> 

have solutions (u, u) G Bt{uo;5]'^ and there is t g]0,T] such that 

(5 2) " ^" ~ > C ( II h - h ||i2(Rd) 

+ /o II (f-f)(s>-) ||L2(]Kd) ds) . 

Inequalities as (5.2) are well-known. In general [7]-[15]-[17], the demonstra- 
tion is achieved in two steps. First detect equilibria where instability arises 
in the discrete spectrum. Then establish that linearized instability implies 
non linear instability. The procedure we adopt below is different. We just 
look at approximate solutions like u^. It follows a more simple proof of (5.2). 

Proof of the Proposition 5.1 . Take I = 2 and N > (8 -|- d). Consider two 
deals of initial data 

iJl{Q,x,e), v\{G,x), i<k<N, 

Ui{0,x,6), ^l{0,x), l<k<N. 
Fix these expressions in the following way 

ul{o,-) = uUo,-), ^\{o,-) = 'pI{o,-) = o, ^Ko,-) = <^i(o,-) = o. 

It implies that 

Ul{t,-)^UHtr), ip\{t,-) = 'pl{t,-) = 0, yte[o,T]. 

Adjust U^{0, ■) and C/|(0, •) so that 

dt {<pI - ifilM •) = - Vifo ■ {Ul - Ui){0, •) ^ . 

Therefore, we are sure to find some t > such that ((^2 ~ '/'2)(*5 ') ^ 0. It 
follows that 

ulit,x,d) = ul{t,x,e + <plit,x)) 
^ ' ' # uKt,x,e) = ul{t,x,e + cpl{t,x)). 

Note u^-*^ and u^^ the approximate solutions built with the profiles {U^}k 
and {U^}k- The associated error terms are f^^ and f^^. 
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Proceed by contradiction. Suppose that the Proposition 5.1 is wrong. Then, 
there is C > and £i G ]0, £o] such that for all e G ]0, ei], we have 

II (uf - uf •) IIl2(R'*) < C ( II - uf )(0, •) ||,.2(R.) 

+ /o ll(f,^'-fb^')(^,-) ||L^(M<i) ds). 
Divide this inequality by y/e. By construction, we have 

£-5 II (uf -uf)(0,-) ||i2(K.)= O(V^), 

II (ff - ff •) IIl2(rc<) = 0(Vi) , Vs G [0,t] . 

£-5 II (uf -uf)(t,-) ||i2(iid) 

= 11 (C/i^-C/f)(t,-,^-V|(«,-) ||l2(mc<) + 0(Vi)- 

It follows that 

lim £-^ II (uf - uf •) ||^2(K<i) = II {Ul - Uf){t, •) |L2(KdxT) = 



e — vO 



which is inconsistent with (5.3). □ 

Remark 5.1.1: In the demonstration presented above, the amplification is 
due to (^2 which is the principal term in the adjusting phase. The presence 
of (p2 becomes efficient in comparison with the other effects when 

This requires to wait a lapse of time bigger than y/e. This delay can be 
reduced by adapting the above procedure to the cases I > 2. A 

Obvious instabilities have an important consequence. To describe the related 
amplifications, it is necessary to introduce new quantities which correspond 
to the phase shifts. In other words, the only way to get —estimates is to 
blow up the state variables. This principle is detailed in [6] in the case of 
compressible Euler equations. 

• Hidden instabilities. Hidden instabilities arc the amplifications which 
are not detected by the monophase description of the section 4. On the other 
hand, they can be revealed by a multiphase analysis. Introduce a second 
phase il)Q{t,x) G WjP such that 

dti>o + (uo • V) V'o = , VV-o A V(/?o ^ 
and disturb the Cauchy data of (5.1) according to 

u^(0, x) = u^(0, x) + ef U{x, £-^ Vo(0, x)) , M :^ N . 
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M 

The small oscillations contained in the perturbation of size e i are not always 
kept under control. They interact with and with themselves. They can 
be organized in such a way to affect the leading oscillation u^. Concretely 
(see [7]), we can adjust U and Vo so that there is a constant C > and 
times ts E]0,Ti;[ going to zero with s such that 

II (u^-Ub)(*'-) \\lHr^)> C VeG]0,eo]- 

The power ef at the time t = is turned into £2 at the time t = t^. 
Such amplifications occur whatever the selection of Z > 2. They imply 
minorations like (5.2). However, the underlying mechanisms are distinct 
from the preceding ones. They are implemented by oscillations which are 
transversal to (^0 and whose wavelengths are OC^)- They are cancelled by 
the addition of the anisotropic viscosity u E^^. 

5.2 Exact solutions 

• Statement of the result. The first information brought by the BKW 

construction is that mean values Uk and oscillations of the profiles Uk 
do not play the same part. This fact is well illustrated by the rules of 
transformation (3.5). It means that we have to distinguish these quantities 
if we want to go further in the analysis. This can be done by involving 
the variables {t,x,9) that is by working at the level of (4.26). To deal with 
{u^ ,p^){t, X, 9) instead of (u*^, p^)(t, x) is usual in non linear geometric optics 
[25]. It allows to mark the terms apt to induce instabilities. 

Select some approximate solution {u^^^ , p^^^ n) ) '^it^ source term fj^ 
given by the Proposition 4.2 and look at 



(5.4) 



f o,£ + (""^ • ST^ciSp u"^ + gtaD^ = v e AvF , 5it)^ = , 
u^{{),x,e) = n|2,^)(0,x,^). 



Theorem 5.1. Fix any integer N > d + 8. There is £n &]0,l] and un > 
such that for alls G]0,eAr] and for all v > the Cauchy problem (5.4) has 
a unique solution {u^,p^) defined on the strip [0,T] x M"' x T. Moreover 



\2,N). 

Proof of the Theorem 5.1 . The system (5.4) amounts to the same thing as 



(5.5) 



dtv!^ + {u^ ■ V)tt^ + div {u^* (g) u^*) + Vp^ = u AxU' 
do,eU^* + [vF ■ gtaDp u^* + e {u^* ■ V) 

+ [(it^* • 0rat)p u^*Y + graD^p^* = u e Au^ 



div = livol = . 
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The equation (5.5) is also equivalent to solve the Cauchy problem 
(5.6) 



' P dtvf + P [(vf ■ V)vf] + P [div (n=* ® u^*)] =u A^u^, 
^ do,eu'* + ^ [{u' ■ mK) u'*] + e ^ [{u'* ■ V) u'] 

associated with the compatible initial data 

u%0, •) = Pu^iO, •) , u'*{0, •) = (0, •) . 

• Blow up. Introduce the new unknown 

:= ( £-T (u- - ul) , {u^* -ul*)), b = (2, N) . 

This transformation agrees with (3.5). The weight e~T in front of (n^ — up 
induces a shift on the indice I. Functions Ui and U*_^ play now the same part 
related to the amplifications. To write the equation on in an abbreviated 
form, we need notations. Quasilinear terms 



Cl,d := P[iul-V)d], 



Cl^d* : 

£•22 d* ■ 



.£* 



+ d* 



I £ 2 u? 



= P [ div (e- 
e-2 [«-V)J], 
-- ^ [ {ul . V)d* ] + e-^ ^1 [ {ul* . graOp d* ] ' 
+ e-^ % [{dtipl + ul-Vipl)ded*] . 



Semilinear terms 

A\^d:= P[{d.V)ul], 
Al^d:= %[{d.Qxal>l){e- 

d* := % [ {d* ■V)ul]+ e-' % [ {d* ■ gxadl) u^* 

Small quadratic terms 

Ql := el P [div {d^d)] + P [div {d* d*)] , 
Ql := ^ [{d-5tadl)d*] + el ^ [{d*-V)d] 
+ [{d* -gtaT,!) d*Y . 

And error terms 

erf := P , erf := s-''' ^ . 

With these conventions, the expression d^ is subjected to 

' P dtdF + Cl^ d^ + d^* + 1 (f _ 
+ £^-1 Q_f + erl = u P A^f , 
^ dtd^* + ^1 d" + >C|2 d^* + d^ + d'* 
+ e'-^ + = Ad^* . 



(5.7) 
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Energy estimates are obtained at the level of (5.7). Below, we just sketch 
the related arguments which are classical. 

• L^— estimates for the linear problem. The linearized equations of 
Euler equations along the approximate solution are obtained by removing 
Qf and from (5.7). It yields a system which, at first sight, involves 
coefficients which are singular in e. In fact, this is not the case. Let us 
explain why. 

This is clear for £21 ^^id Af-^^. 

Since u^* = C){£~), this is also true for and A21. 

The contributions which in have in factor give no trouble since 

dtv^^ + ul ■ V<^f = 0{e^) = 0{e''+^) , < • Vipl = vl* = 0{e'+J) . 
Now, look at A22- Recall that (F* = ^^d^* which means that 

£-1 d^* • V(^f = -divd^*. 
Therefore 

£-1 ^ [{d'*-5xaT)l)ul*]* = T^{t,x,V)d'*, 

where is some differential operator of order 1 with bounded coefficients. 

Observe that these manipulations and the blow up procedure induce a loss 
of hypcrbolicity. When u = 0, this is the source of hidden instabilities. 
When > i^jv > with z/^v large enough, this can be compensated by the 
viscosity. This is the key to L^— estimates. 

• The non linear problem and higher order estimates. Let a be the 

smaller integer such that a > If the life span of the exact solution 
is finite, we must have 

lim II u'^it,-) \\h''= +00. 

t >Te 

Thus, the Theorem 5.1 is a consequence of the following maj oration 

sup { II u^it, •) \\h- ; i G [0, min {T^,T)] } < C < 00 . 
Consider the set 

■= { <^0,e , ■■■ , f d,£ ,90}. 
Extract the operators 

:= Zi o ■ ■ ■ o , Zj E Z^ , k < a . 
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It suffices to show that 



max sup { II u%t, •) 11^2 ; t e [0, min (T^, T)] } < C < oo . 

Pick some with k < a. Apply Z^ on the left of (5.7). Use the point i) of 
Lemma 4.2 to pass through Then, observe that the commutator of two 
vector fields in Z'^ is a linear combination of elements of Z^ with coefficients 
in C°°. Thus, we get an equation on Z^ d^* . 

The linear part is managed as in the preceding paragraph. Take i = 1. The 
contributions due to Qf and Ql are controled by way of the a priori estimate 
and the viscosity. The condition on N is to make sure that 
N , 3 ^ ^ n 

Thereby, the contributions brought by the error terms erf and erf remain 
bounded in the procedure. 



References 

[1] C. Bardos, What use for the mathematical theory of the Navier- 
Stokes equations, Mathematical fluid mechanics, 1-25, Adv. Math. Fluid 
Mech., Birkhauser, Basel (2001). 

[2] J. T. Beale, T. Kato, A. Majda, Remarks on the breakdown of smooth 
solutions for the 3-D Euler equations, Commun. Math. Phys., P. 61-66 
(1994). 

[3] T. Chacon Rebollo, Oscillations due to the transport of microstructures, 
SIAM J. Appl. Math., Vol 48, No 5, Octobre 1988. 

[4] J.-Y. Chemin, Perfect incompressible fluids. Translated from the 1995 
French original by Isabelle Gallagher and Dragos Iftimie. Oxford Lec- 
ture Series in Mathematics and its Applications, 14. The Clarendon 
Press, Oxford University Press, New York, 1998. x-|-187 pp. 

[5] J-Y. Chemin, B. Desjardins, I. Gallagher, E. Grenier, Fluids with 
anisotropic viscosity, Mathematical Modelling and Numerical Analy- 
sis, M2AN, Vol. 34, 2, 315-335 (2000). 

[6] C. Cheverry, Propagation of Oscillations in Real Vanishing Viscosity 
Limit, To appear in Commun. Math. Phys. 



40 



[7] C. Cheverry, O. Gues, G. Metivier, Oscillations fortes sur un champ 
lineairement degenere, Annales Scientifiques de PENS, 5 (2003). 



[8] C. Cheverry, O. Gues, G. Metivier, Large amplitude high frequency 
waves for quasilinear hyperbolic systems, submitted in Advances in Dif- 
ferential equations. 

[9] P. Constantin, C. Fefferman, Direction of vorticity and the problem of 
global regularity for the Navier-Stokes equations, Indiana Univ. Math. 
Journal, 42, No. 3, 775-789 (1993). 

[10] A. Corli, O. Gues, Stratified solutions for systems of conservation laws. 
Trans. Amer. Math. Soc. 353, no. 6, 2459-2486 (2001). 

[11] R.-J. DiPerna, A.-J. Majda, Oscillations and Concentrations in Weak 
Solutions of the Incompressible Fluid Equations, Commun. Math. Phys. 
108, 667-689 (1987). 

[12] W. E, Homogenization of linear and non linear transport equations. 
Comm. Pure Appl. Math. 45, no. 3, 301-326 (1992). 

[13] W. E, Propagation of oscillations in the solutions ofl — D compressible 
fluid equations. Com. P. D. E. 17 (3/4), 347-370 (1992). 

[14] C. Foias, P.-P. Manley, R. Rosa, R. Temam, Cascade of energy in tur- 
bulent flows, C. R. Acad. Sci. Paris, t. 332, Serie I, p. 509-514 (2001). 

[15] S. Priedlandler, W. Strauss, M. Vishik, Nonlinear instability in an ideal 
fluid, Ann. Inst. Henri Poincare, Vol. 14, no. 2 (1997), 187-209. 

[16] P. Gerard, Microlocal defect measures. Comm. in Partial Differential 
Equations, 16, 1761-1794 (1991). 

[17] E, Grenier, On the nonlinear instability of Euler and Prandtl equations. 
Comm. Pure Appl. Math. 53 (2000), no. 9, 1067-1091. 

[18] O. Gues, Developpement asymptotique de solutions exactes de systemes 
hyperboliques quasilineaires. Asymptotic Anal. 6, no. 3, 241-269 (1993). 

[19] O. Gues, Ondes multidimensionnelles e—stratifiees et oscillations, Duke 
Math. J. 68, no. 3, 401-446 (1992). 

[20] J-L Joly, G. Metivier, J. Ranch, Transparent non linear geometric optics 
and Maxwell-Bloch equations. Journal of Differential Equations, 166, 
175-250 (2000). 



41 



[21] D. Lannes, Secular growth for symmetric hyperbolic systems, Journal of 
Differential Equations, 90, 466-503 (2003). 

[22] M. Lesieur, Turbulence in fluids, Fluid Mechanics and its Applications, 
Kluwer Academic Publishers Group, Dordrecht, 515 pp. ISBN: 0-7923- 
4415-4 (1997). 

[23] P.-L. Lions, Mathematical Topics in fluid Mechanics, Volume 1, Incom- 
pressible Models, Oxford Science Publications. 

[24] D. W. Mclaughlin, G. C. Papanicolaou, O. R. Pironneau, Convection 
of micro structure and related problems, SIAM J. Appl. Math., Vol 45, 
No 5, Octobre 1985. 

[25] S. Schochct, Fast singular limits of hyperbolic PDEs, J. Dif. Eq. 114, 
476-512 (1994). 

[26] D. Serre, Oscillations nonlineaires de haute frequence ; dim > 2, In 
Marino A. and Murthy M. K. V., editors. Nonlinear variational prob- 
lems and partial differential equations, volume 320 of Pitman Res. notes 
in Math., pages 245-294, London, 1995. Longman. 

[27] T. Sideris, Formation of Singularities in Compressible Fluids, Gommun. 
in Math. Phys. 101, 475-485 (1985). 



42 



